Kinetic equation for a dense soliton gas 
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We propose a general method to derive kinetic equations for dense soliton gases in physical 
systems described by integrable nonlinear wave equations. The kinetic equation describes evolution 
of the spectral distribution function of solitons due to soliton-soliton collisions. Owing to complete 
integrability of the soliton equations, only pairwise soliton interactions contribute to the solution 
and the evolution reduces to a transport of the eigenvalues of the associated spectral problem with 
the corresponding soliton velocities modified by the collisions. The proposed general procedure of 
the derivation of the kinetic equation is illustrated by the examples of the Korteweg - de Vries 
(KdV) and nonlinear Schrodinger (NLS) equations. As a simple physical example we construct an 
explicit solution for the case of interaction of two cold NLS soliton gases. 

PACS numbers: 05.45.Yv 



The concept of soliton plays a fundamental role in non- 
linear physics due to its main property — preservation of 
its parameters during interactions with other solitons in 
the case of the physically and mathematically important 
class of so-called integrable equations. Such a particle- 
like behavior has led to a huge number of investigations 
of soliton dynamics in various physical systems (see, e.g. 
0) as well as to thorough study of mathematical prop- 
erties of integrable (or soliton) equations (see, e.g. |2|) 
initiated in the pioneering paper |3| where the famous in- 
verse scattering transform (1ST) method has been formu- 
lated. In this method, each soliton is parameterized by 
an eigenvalue A of the linear spectral problem associated 
with the nonlinear wave equation under consideration. 
For example, the KdV equation 



ut + 6uu x + u xx 







is associated with the linear Schrodinger equation 



ip xx + u(x,t)%jj 



(i) 



(2) 



so that evolution of the potential u(x,t) according to 
does not change the spectrum A, and, whence, the prop- 
erties of solitons do not change either. The 1ST method 
gives a full explanation of finite-number solitons dynam- 
ics and provides the basis for the description of many 
nonlinear physics phenomena. 

We encounter a different physical situation when we 
have to deal with a dense lattice of solitons. When the 
solitons in the lattice are correlated they may form a 
modulated nonlinear periodic wave. The spectral prob- 
lem (J2J) for a general periodic in x "potential" u(x, t) leads 
to the Bloch band structure of the spectrum A. The 
periodic "soliton lattices" are distinguished by a finite 
number of gaps in the spectrum (see 0) and the edges 
Aj of the spectral gaps become convenient parameters in 
terms of which the major physical characteristics of the 
wave such as wavelength, frequency, amplitude etc are 
expressed. In a weakly modulated wave the parameters 



Ai become slow functions of space and time coordinates, 
and their evolution is governed by the Whitham equa- 
tions Methods of derivation and integration of the 
Whitham equations are now well developed and this the- 
ory provides an adequate description of such important 
phenomenon as formation of dispersive shocks (or undu- 
lar bores) in various physical systems from water surface 
to space plasma and Bose-Einstein condensate. 

Yet another principally different class of problems 
arises when solitons form a disordered finite-density en- 
semble (a soliton gas) rather than well-ordered modu- 
lated soliton lattice. The relevant physical conditions 
can be realized for instance when large number of soli- 
tons are generated either by random large-scale initial 
distributions 0] or by a stochastic external forcing (such 
as irregular topography in the internal water wave dy- 
namics), or when they are injected into a ring resonator 
Here it is necessary to introduce an appropriate 
kinetic description of the soliton gas. We introduce a 
distribution function f(x,t; A) as the number of solitons 
with the spectral parameter A in the interval (A, A + dX) 
and in the space interval (x, x + dx) at the moment t. 
Now, if the soliton dynamics is governed by an exactly 
integrable equation we arrive at the problem of describ- 
ing the isospectral evolution of the distribution function 
f(x,t;X) with time. Since the spectrum A is preserved, 
the evolution of f(x,t;X) must be governed by the con- 
servation equation 



A + (*/)* = o ) 



(3) 



which means that the eigenvalues A are transferred along 
the a;-axis with some mean velocity s(x, t; X) depending 
on the distribution function / but there is no exchange 
of A's between different A-intervals. Hence, the problem 
is reduced to finding the velocity s of a soliton gas as 
a function of A, x and t. This problem was posed by 
Zakharov || as early as in 1971 and he solved it for the 
case of rarefied gas of KdV solitons. Only recently it 
was generalized in Ref. |9j] to the case of a dense gas of 



2 



KdV solitons where it was found that velocity s(x, f; A) 
is determined from the integral equation 



s (A) = 4A 2 + i J 



In 



A + /t 



f(n)[s(\) - s(ti)}dfi, (4) 



so that Eqs. © and J3J give a closed self-consistent ki- 
netic description of a soliton gas of an arbitrary density. 
In the limit of rarefied gas, i.e. for J f(p)dfj, <C Ao, where 
Ao is a characteristic value of the spectral parameter, the 
second term in Eq. (0} becomes a small correction to the 
speeds 4A 2 of non-interacting solitons, and the substitu- 
tion of s(r/) = 4t? 2 , (rj = A, /i), reproduces the Zakharov 
kinetic equation |s|. 

Although a mathematically rigorous derivation of 
Eq. given in Ref. (which is based on a certain 
singular limit of the Whitham equations) is quite tech- 
nical, the final result is physically very natural and sug- 
gestive. Indeed, Eq. Q implies that only two-soliton 
collisions have to be taken into account (which agrees 
with the properties of multi-soliton solutions of the KdV 
equation); thus in a collision of a A-soliton (i.e. with the 
spectral parameter A) with a /i-soliton, the coordinate of 
the A-soliton is shifted by the distance 

A > H, 



A 



for 



(and there is a similar expression for A < //), and num- 
ber of collisions per second is proportional to the relative 
mean velocity [s(fi) — s(A)] of these two types of solitons 
multiplied by the density of /x-solitons. Then, after inte- 
gration over the distribution function /(/z) of /z-solitons, 
we arrive at equation Q for the speed of A-solitons mod- 
ified by their collisions with the other /i-solitons. It is 
supposed that the number of the soliton collisions over 
large distance is large enough and, hence, the mean ve- 
locity s is a well-defined variable. As a matter of fact, 
the typical x, i-scales in the kinetic equation are much 
larger than in the original equation QJ. Thus, one can 
see that the above simple reasoning provides an indepen- 
dent derivation of Eq. (@J for the KdV equation case and, 
obviously, it can be directly applied to other integrable 
equations. In this Letter we use this method to derive 
the kinetic equation for a finite-density gas of bright NLS 
solitons and infer some its consequences. 

As is known, each soliton solution of the focusing NLS 
equation 



iu t 



2\u\ z u = 



(5) 



is characterized by a complex eigenvalue 

A = a + 17, — oo < a < oo, < 7 < 00, (6) 

of the Zakharov-Shabat spectral problem and is given by 
(see, e.g. 0) 

, o- exp[-2iax - 4i(a 2 - j 2 )t - i(j} ] 

u(x,t) = 2i r y , . . - , (7) 

y ' ' ' cosh[2 7 (a; + 4at - x )} 



that is 7 = 3(A) determines the amplitude of the soliton 
and a = 3?(A) its velocity v — —4a. The multi-soliton so- 
lution shows that interaction of solitons reduces to only 
two-soliton elastic collisions and effects of multi-soliton 
collisions are absent. In accordance with the above argu- 
mentation, we consider a gas of solitons characterized by 
a continuous distribution function of eigenvalues A. In 
other words, f(x,t]a,^)dad'ydx is the number of eigen- 
values in the element dad-y of the complex plane A in the 
space interval dx much greater than both typical soliton 
width ~ 1/7 and average distance between solitons. This 
distribution function evolves due to motion of solitons 
and as a consequence of preservation of the spectrum it 
satisfies again the continuity equation © which means 
conservation of the density of eigenvalues in a given ele- 
ment of A-plane. In this equation s denotes mean velocity 
of A-solitons which should be evaluated with the account 
of soliton collisions. Without interaction of solitons it 
would be equal to s(a, 7) = —4a. However, collisions of 
A-soliton (A = a + 17) with other /i-solitons (// = £ + iif) 
modify it in the following way. Each collision of faster A- 
soliton with slower /x-soliton shifts the A-soliton forward 
to the distance (see, e.g. Q) 
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In 



A 



A — /j, 



for s(a,7) > s(£,rj), 



and the number of such collisions in the time interval 
dt is equal to the product of the density /(£, rj)d^di] and 
the distance overcame by faster A-solitons compared with 
slower ^-solitons, [s(a,7) — s(^,rj)]dt. Thus, such colli- 
sions increase a path covered by A-soliton compared with 
—Aadt. In a similar way, one can calculate the nega- 
tive shift due to collisions with faster /i-solitons. The 
total shift is obtained by integration over d^di]. Equat- 
ing paths s(a,j)dt and —Aadt+ "total shift", we obtain 
an integral equation for the self-consistent definition of 
the soliton gas velocities, 



s(a,7) = -4a + — 
2 7 



00 J 



In 



A — /! 



A 



(8) 



x[s(a,7) - s(£,r))]d(dr). 



Obviously, this derivation of Eq. © is correct as long as 
the resulting velocity s(a, 7) is finite. Under this reserva- 
tion, equations 10 and JSJ provide the basis for a consid- 
eration of kinetic behavior of dense gas of NLS solitons. 

As a simple application of the kinetic equations, let 
us consider evolution of two beams of solitons, when the 
spectral distribution gas consists of two monochromatic 
("cold") parts: 



f(x,t;£,v) = Pi(%,t)d(-£ - a,?/ - 7) 
+ p 2 (x,t)5(£ - a,rj - j), 



(9) 



where pi(x, t) corresponds to fast (or moving to the right 
in the reference system associated with the group velocity 
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of the carrier wave) solitons (Ai = —a + £7) and P2{x, t) 
to slow (or moving to the left) solitons (A2 = a + i-f). All 
solitons have the same amplitudes. Of course, the ideal- 
ized delta-functional approximation for the distribution 
function in the kinetic equation means that the exact soli- 
ton eigenvalues in the original spectral problem for the 
NLS equation are distributed for i-th component within 
a narrow vicinity of the dominant value A = Ai and are 
actually different, which precludes formation of bound 
states. The soliton positions in such a "monochromatic" 
gas are statistically independent which leads to Poisson 
distribution for the number of solitons in a unit space 
interval with the mean density pi [ldj . It is also clear 
that one can neglect interactions between the solitons be- 
longing to the same beam compared with the cross-beam 
interactions. 

Substituting Jj|J| into @, (JSJ) we obtain a "two-beam" 
reduction of the kinetic equation, 



dpi 
dt 



djsipi) 
dx 



= ^ 
' dt 



9{s 2 p2) 
dx 



= 0, 



(10) 



where the velocities of the beams s\ = s(— a, 7) and S2 
s(a,"f) are determined from the equations 



si 



4a + np 2 -{si-s 2 ), s 2 — -4a + np 1 ■ (s 2 — Si), (11) 



and the interaction parameter 
1 



2 7 ' 



In 1 



(12) 



is always positive. Resolving Hll|) we get expressions for 
velocities in terms of p\^ 2 (x 1 t): 



si = 4a 



1 - Kjfli - P2 ) 

1 - k(pi + p 2 ) ' 



s 2 = —4a 



1 + k(px - p 2 ) 
1 - k(pi + p 2 )' 



(13) 



According to the formulated above condition of applica- 
bility of the kinetic description, the densities must satisfy 
the inequality 



n(pi + p 2 ) < 1. 
Substitution of the inverse expressions for pi, 2 , 
s 2 + 4a s% — 4a 



Pi = 



k(s 2 - si) ' 2 k(si - s 2 ) 
into the system (|10[) reduces it to the form 



dsi 
~dt 



s 2 



dsi 
dx 



0. 



ds 2 
dt 



ds 2 
si— = 0, 
dx 



(14) 



(15) 



(16) 



which is known as the Riemann invariant form of a Chap- 
lygin gas equations which, besides well-known original 
application to compressible gas dynamics, has recently 
found a number of other applications. If solution of the 
system (|16|) is known, then the densities p\ , p 2 are deter- 
mined in terms of s±, s 2 by the formulas (|15|) . 



Although the system l|16(l admits the general solution 
(see, e.g. QJJ), we shall confine ourselves to a physically 
natural problem of the collision (mixing) of slow and fast 
soliton gases when two gases are separated at the initial 
moment so that 



pi(x, 0) = p w {x), p 2 {x,0) = 



for x < , (17) 



p 1 (x,0) = 0, p 2 {x,0) = p 20 (x) for x > (18) 

where pio(x) > and p 2 o(x) > are given functions. 
Corresponding initial data for the equivalent system 1|16|) 
follow from (|13fl : 



si(x, 0) = 4a, s 2 (x, 0) 



-4a 



1 + Kp 10 (x) 
1 - Kp w (x) ' 



(19) 



for x < 0, and 



si(x, 0) = 4a 



1 + Kp 20 (x) 
1 - Kp 20 (x) ' 



s 2 (x,0) = -4a, (20) 



for x > 0. To further simplify the problem, we con- 
sider the case when both gases are homogeneous, i.e., 
Pw{x) = pio, P2o{x) = P20, where p w and p 20 are con- 
stant. Since in this case neither the initial conditions, 
nor the evolution equations l|ltj|) contain any parame- 
ter with dimension of length, the solution must be self- 
similar and depend on the variable x/t alone. However, 
as can be easily seen, the system (tTfill does not possess 
non-constant similarity solutions. On the other hand, the 
solution consisting of two constants s\ , s 2 cannot sat- 
isfy the discontinuous initial conditions i|19|l . (|20fl . As 
in shock wave theory this can be remedied by in- 
troducing admissible discontinuities in the solution. The 
discontinuous "weak" solutions are allowed here owing 
to the presence of the conservation laws (|10|) , l|13|l . As a 
result, the sought solution has the form of three constant 
states {/i,/ 2 } : {pio, 0}, {pi c , p 2c } 7 {0, p 20 } separated by 
two jump discontinuities. 

We denote the velocities of the discontinuities as 



c , > c 
(i) x < c~t: 



Then the solution is: 



Pl = p x = pio , P2 = P 2 = , 

which implies by (|13|) . 

1 + npia 



Si = — 4a , s 2 = s 2 



-4a- 



1 - K/OlO 



(ii) x > c+t: 



Pl=Pt=°i P2=P 2 =P2Q, 



which implies 



si = sT = 4a , s 2 — s, = —4a. 

1 - Hp 2 



(21) 



(22) 



(23) 



(24) 
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The values and can be viewed as velocities of trial 
solitons of one component moving through the homoge- 
neous gas of solitons of another component. One can see 
that there are critical values for the densities p%o = n^ 1 
and P20 = n^ 1 yielding the infinite speeds for trial soli- 
tons. In accordance with the restriction described above, 
we assume npio < 1, KP20 < 1. 
(hi) c~t < x < c + t: 

Let pi = pic, P2 = P2c in this interaction zone. Then 
the four unknown values p\ c , p2c, c + , c~ are found from 
the jump conditions consistent with the physical conser- 
vation laws iJTUj l (see, for instance, ||): 



[Pic 
[P2c 



P2 



-c + [pic - pi) 

-C + [p 2 c - Pt) 



[picSlc 
[p2cS2c 



[PicSlc 
[p2cS2c 



Pljl_ 
P2 S 2 



= 0, 
= 0, 



0, 



pf4 
pt4\ = 



(25) 



(26) 



Here 



sic = 4a 



1 - k(pi 



P2c 



1 - n(pic + p 2c ) ' 



*2c 



-4a 



1 + n(pi c - p 2c ) 

1 - K(pi c + p 2c )' 

(27) 

In view of J23 it readily follows from and lf^|) 
that c~ = S2c and c + = s± c . Then p\ c , p2 C are found 
from the same equations with the account of (|2"T|) . i(2"3"|) 

as 



Pic 



Plo(l - KP20) 
1 - K 2 p w p20 ' 



P2c 



P20(l - KPw) 



1 



K^PwP20 



(28) 



Correspondingly, substitution of (|28|l into l|27(l yields the 
velocities of expansion of the interaction region 



S2c 



-4a- 



1 



K-PlO 



1 - Kp 10 



Sic 



Act - 



1 



1 - ' 



(29) 

and this completes the solution. As one should expect, 
these velocities coincide with velocities and s\ of the 
trial solitons (see Eqs. Il22[l and J2U)- Thus, we have 
found boundaries x — c f of the interaction region and 
densities pi c , /?2c of the two components of soliton gas in 
this region. In particular, the obtained solution in view 
of (I14|) implies that the following inequality is satisfied: 



F(pi c ,P2c 



1 



2^:piop2Q 
P10+P02 



Pic + P2c 
PlQ + P20 1 - H 2 p 10 p 2 



< 1. 



(30) 



It means that the total density in the interaction (mixing) 
zone is always less than the sum of densities of individual 
separated components. 

The process of the collision of two soliton gases is il- 
lustrated in Fig. 1. Due to the interaction of solitons 
with each other, the overlap region spreads out faster 
than it would without taking into account the phase- 
shifts caused by two-soliton collisions, and the kinetic 





Plc+P 2 c 






P10 


Pic 




P 20 






P2c 









c"t 



c+t 



FIG. 1: Dependence of solitons densities on space coordinate 
x at some moment t; pio and p2o are the densities of solitons 
in the beams propagating to the left and to the right, corre- 
spondingly; pic and p2c are the densities of the same beams 
in the overlap region. 



equation allows one to give quantitative description of 
this effect. 

In conclusion, we have obtained the kinetic equation 
describing the evolution of the spectral distribution func- 
tion of a dense gas of uncorrelated NLS solitons. The 
proposed procedure of the derivation can be generalized 
to whole AKNS hierarchy and to other integrable hierar- 
chies. A process of interaction of two "cold" NLS soliton 
gases is studied in detail. 
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